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ABSTRACT 

Aims. We model the synthesis of molecules and dust in the inner wind of the oxygen-rich Mira-type star IK Tau by considering the 
effects of periodic shocks induced by the stellar pulsation on the gas and by following the non-equilibrium chemistry in the shocked 
gas layers between 1 R* and 10 R*. We consider a very complete set of molecules and dust clusters, and combine the nucleation phase 
of dust formation with the condensation of these clusters into dust grains. We also test the impact of increasing the local gas density. 
Our derived molecular abundances and dust properties are compared to the most recent observational data. 

Methods. A semi-analytical formalism based on parameterised fluid equations is used to describe the gas density, velocity, and 
temperature in the inner wind. The chemistry is described by using a chemical kinetic network of reactions and the condensation 
mechanism is described by a Brownian formalism. A set of stiff, ordinary, coupled differential equations is solved, and molecular 
abundances, dust cluster abundances, grain size distributions and dust masses are derived. 

Results. The shocks drive an active non-equilibrium chemistry in the dust formation zone of IK Tau where the collision destruction 
of CO in the post-shock gas triggers the formation of C-bearing species such as HCN and CS. Most of the modelled molecular 
abundances agree well with the latest values derived from Herschel data, except for SO2 and NH3, whose formation may not occur in 
the inner wind. Clusters of alumina, AI2O3, are produced within 2 R* and lead to a population of alumina grains close to the stellar 
surface. Clusters of silicates (Mg2Si04) form at larger radii (r > 3 R*), where their nucleation is triggered by the formation of HSiO 
and H2SiO. They efficiently condense and reach their final grain size distribution between ~ 6 R* and 8 R* with a major population 
of medium size grains peaking at ~ 200 A. This two dust-shell configuration agrees with recent interferometric observations. The 
derived dust-to-gas mass ratio for IK Tau is in the range 1 - 6 x 10“^ and agrees with values derived from observations of 0-rich 
Mira-type stars. 

Conclusions. Our results confirm the importance of periodic shocks in chemically shaping the inner wind of AGB stars and providing 
gas conditions conducive to the efficient synthesis of molecules and dust by non-equilibrium processes. They indicate that the wind 
acceleration will possibly develop in the radius range 4 - 8 R* in IK Tau. 

Key words. Stars: AGB and post-AGB; (Stars:) circumstellar matter, Astrochemistry; Stars: Individual: IK Tau; ISM: dust, extinction; 
Molecular processes 


1. Introduction 

Low mass stars with masses between 1 Mq and 4 Mq ascend the 
asymptotic giant branch (AGB) in the late stages of their evo¬ 
lution where they experience thermal pulses, dust formation in 
their cool extended atmosphere and mass loss. AGB stars are the 
prevalent dust providers to galaxies through their strong winds 
that replenish the surrounding medium with molecules and dust 
grains (IGehrz 19891) . Depending on the number of dredge-up 
events experienced by the star and its evolutionary stage on the 
AGB, the red giant forms either oxygen-rich or carbon-rich solid 
grains reflecting the carbon-to-oxygen (hereafter C/O) ratio of 
the stellar photosphere. The astronomical evidence for the syn¬ 
thesis of dust in AGBs is brought by the identification of typ¬ 
ical dust features in the spectral energy distribution (SED) of 
AGB stars, including the 9.7 pm and 18 pm features of sili¬ 
cates in oxygen-rich AGBs and the 11.3 pm feature ascribed to 
silicon carbide grains in carbon stars. Recent surveys with the 
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Herschel space telescope of the Magellanic Clouds have ques¬ 
tioned the primacy of AGBs as galactic dust providers at lower 
metallicity than that of the Milky Way jMatsuura et al. 20091 
Boyer et al. 2012| l. However, owing to their large number in stel¬ 
lar populations, AGB stars remain important dust producers in 
the local and high redshift universe. 

Despite their crucial role in the dust budgets of galax¬ 
ies, the mechanisms responsible for the production of dust in 
AGB winds are still largely unknown. Erom interferometric 
data, we know that the dust formation region is located within 
~ 10 R* from the stellar photosphere dPanchi et al. 1994l l. 
Periodic shocks resulting from the stellar pulsation per¬ 
vade the region, lift dense and warm gas layers to larger 
radii, and facilitate the formation of solids and the accel¬ 
eration of the stellar wind by radiation pressure on dust 
grains because of their high dust opacity (IBowen 1988t . 
These dense layers are also conducive to the formation of 
molecules, and many gas-phase species have been identified 
close to the stellar surface in both oxygen- and carbon-rich 
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AGBs by various observing facilities, including most recently 
Herschel Ie.g.. lCernicharo et al. 20081 Groenewegen et al. 2011 


ISchoi'er et al, 201 fl lJusttanont et al. 20121 . Furthermore, these 
layers are characterised by non-equilibrium chemical processes 
induced by the periodic passage of the shocks that leads to 
the formation of molecules and dust grains. For molecules, 
the chemistry in the photosphere is controlled by the C/O ra¬ 
tio and the formation of carbon monoxide, CO. This molecule 
traps all oxygen (carbon) atoms in a carbon-rich (oxygen- 
rich) gas. The partial dissociation of CO by the shocks re¬ 
leases oxygen atoms in carbon-rich environments and car¬ 
bon atoms in oxygen-rich environments, leading to a dual 
chemistry that forms O-bearing species like water, H 2 O, in 
the wind of carbon stars, and C-bearing species in the wind 
of O-rich AGBs ( [Willacy & Cherchneff 199^IDuari et al.~7999l 
Cherchneff 2006, 2011a, 2012). The widespread occurrence 
of warm water in several carbon stars has been highlighted 
with Herschel (INeufeld et al. 20lT1) whereas hydrogen cyanide, 
HCN, was detected close to the stellar surface in O-rich Miras, 
again with Herschel (lJusttanont et al. 20121 . These detections 
confirm the crucial role of shock-induced chemistry in the for¬ 
mation of these unexpected molecules in AGB stellar winds, al¬ 
beit other processes have been proposed that may also contribute 
in breaking the chemical dichotomy ( Aguhdez et al. 20T0| . 

The SEDs of oxygen-rich AGB stars show spectral bands 
typical of alumina (AI 2 O 3 ) in the range 9-15 pm, of sili¬ 
cates around 9.7 pm and 18 pm, or a mix of both. Furthermore, 
stars that ascend the AGB show trends in dust production 
dLittle-Marenin & Little 1990] ). The less evolved, low mass-loss 
rate stars show the presence of alumina (AI 2 O 3 ), while more 
evolved stars characterised by high mass-loss rates show sili¬ 
cate features. Observations point to the synthesis of large grains 
{a > 0.3 pm) close to the star in a few O-rich objects that in¬ 
clude semi-regular variable stars (INorris et al. 2012t . while the 
presence of alumina and/or silicates close to the star (r < 2 R*) 
and at larger radius {r ~ 5 R*), respectively, have been inferred 
to explain recent interferometric data (IWittkowski et al. 20071 
ISacuto et al. 20131 IKarovicova et al. 20131) . Because the alu¬ 
mina and silicate opacities are low in the wavelength range 
where most of the stellar flux is emitted, scattering on large 
grains {a > 0.1 pm) was proposed as a viable mechanism to 
accelerate the wind in O-rich AGBs at small radii dHofner 20081 
IBladh et al. M 2 I) . 

Despite the astronomical evidence of alumina and sili¬ 
cates in O-rich AGBs, the processes controlling the forma¬ 
tion of dust in these stars are not well understood. Previous 
theoretical studies on dust production in O-rich AGB stars 
have relied on homogeneous or heterogeneous nucleation the¬ 
ory. They usually assume that chemical equilibrium prevails 
in the dust formation zone, and derive the gas phase compo¬ 
sition from equilibrium conditions. They further assume that 
the growth of dust grains proceeds through surface deposition 
of atoms and molecules, and the final size of grains is con¬ 
trolled by accretion and evaporation under equilibrium con- 
ditions ([Gail & Sedlmayr 1999 Ferraroti & Gail 2002, 2006, 
IWoitke 20061 iHofner 2008r De I’Agli et al. 2014| l. Some inves¬ 
tigations have pointed to the fact that silicates and oxides such 
as alumina do not have gas-phase homogeneous monomers out 
of which condensation of grains proceeds, and that heteroge¬ 
neous condensation may occur through the deposition of chem¬ 
ical species on the surface of pre-existing grain seeds of dif¬ 
ferent chemical type (e.g., alumina, AI 2 O 3 , or titanium ox¬ 
ide, Ti 02 ) (Gail & Sedlmayr 1998 jJeong et al. 2003 Ferraroti 
& Gail 2002, 2006). In some models that aim to explain wind 


acceleration and reproduce the photometric properties of O-rich 
AGBs, a pre-requisite to dust growth by surface deposition is the 
assumption that a population of small dust seeds exists (Bladh 
et al. 2012, 2013, 2015). Neither the formation of these dust 
grains nor the formation of molecules from the gas phase by 
non-equilibrium processes is treated in such approaches. A thor¬ 
ough review of these various topics and mechanisms involved 
in homogeneous and heterogenous nucleation at chemical equi¬ 
librium is provided by Gail & Sedlmayr (2013). Recently, the 
role of titanate clusters as seed nuclei in AGB winds was investi¬ 
gated (IPlane 20131|Goumans & Bromley 2013| l, while Goumans 
& Bromley (2012) studied the formation pathways to silicate 
dust from the gas phase. In this latter study, the growth of en- 
statite, MgSi 03 , and forsterite, Mg 2 Si 04 , clusters follows the 
combination of oxygen addition and magnesium inclusion steps. 
Depending on the gas conditions, these chemical routes become 
exergonic and lead to the efficient synthesis of enstatite and 
forsterite dimers. Whatever the chemical processes involved in 
the synthesis of dust are, it is likely that the dense, shocked lay¬ 
ers above the photosphere of AGB stars provide the ideal gas 
conditions for the synthesis of molecules, dust clusters and the 
growth of dust grains by non-equilibrium processes (Cherchneff 
2012, 2013). 

Building on these findings and on a previous study by Duari 
et al. (1999), we present in this paper a complete chemical de¬ 
scription of the inner wind of the oxygen-rich AGB star, IK Tauri 
(or NML Tauri). This exhaustive model describes the formation 
of molecules and dust grains in the inner wind of an oxygen- 
rich AGB star and is based on a non-equilibrium, chemical ki¬ 
netic approach. We consider the impact of the periodic shocks on 
the gas layers above the photosphere and study the gas and the 
solid phase of the extended atmosphere, which corresponds to 
the dust formation and wind acceleration zone between 1 and 
10 R*. We model the formation of various chemical species, 
some of which have being recently detected in the inner wind 
of IK Tau with Herschel, and simultaneously model the forma¬ 
tion of oxygen-rich clusters, including silicates and alumina, out 
of the gas phase. We then follow the condensation of these clus¬ 
ters to form dust grains and derive grain size distributions. This 
study aims to derive trends on the molecular component and the 
formation of dust in the inner wind of IK Tau, but the results 
may equally apply to other O-rich stars in the Mira phase on the 
thermally-pulsing AGB. 

The paper is structured as follows; we briefly present the 
physical model used to describe the shocked inner wind of IK 
Tau in ^ The chemistry of molecule and cluster formation, and 
the condensation formalism are described in 92. 2| Results on 
molecules, clusters and dust are presented in ^ and finally, a 
discussion follows in ^ 

2. Model of the inner wind 

The photosphere is assumed to be at thermodynamic equilibrium 
(TE) and is characterised by one gas temperature and pressure at 
which molecular abundances are estimated. We then shock the 
gas layer at the photosphere, which gradually moves to larger 
radii where it is again shocked, as described in § [2.1 j We as¬ 
sume spherical symmetry of the dust formation zone, and that 
one shock propagates in the inner wind at each pulsation pe¬ 
riod. In the chemical description of the gas phase, we consider 
all molecules detected close to the stellar surface in O-rich Miras 
stars, and the formation of dust clusters that triggers the conden¬ 
sation of dust grains. We consider two chemical types of dust, 
silicates and alumina, and the molecular clusters that partake in 
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Fig. 1. Schematic view of the dust formation zone in O-rich AGB stars, which includes typical physical parameters and the prevalent 
chemical processes related to dust production. Molecules present under TE in the photosphere are shown. 


Table 1. Stellar and inner wind parameters for IK Tau. The wind 


parameters are explained in § 2.1 


Star 

M*,zams 

1.5 Mo 

1 

R* 

2.5 xlO‘3 cm 

2 

C/O 

0.75 

3 

P 

470 days 

4 

Te// 

2200 K 

2 

Wind 

r* 

IR* 

3 

V, 

32 km s“' 

3 

M 

8 X lO-* Mo yr-‘ 

2 

''^esc 

17.7kms-‘ 

2 

7 

0.98 


a 

-0.6 

3 

nga.,(r.5) 

3.62x10'^ cm-3 

3 

Ts<,.Tr.,) 

2200 K 

- 


References. 1) Fox & Wood (1982) ; 2) Decin et al 2010; 3) Duari et 
al. (1999); 4) Hale et al. (1997) 


the formation of small condensation nuclei. We condense the 
dust seeds as soon as they form, by assuming they coalesce and 
coagulate to form dust grains. The studied region spans 1 R* to 
10 R*, and is illustrated in Figure where the gas parameters 
and the main chemical and physical processes are indicated. 

2.1. Physical model 

To describe the physical conditions experienced by the up¬ 
per atmosphere of IK Tau, lifted by the recurrent passage 
of shocks induced by the stellar pulsation, we use the semi- 
analytical model of the inner wind of AGB stars developed by 


Cherchneff et al. (1992) and Cherchneff (1996), and used by 
Willacy & Cherchneff (1998), Duari et al. (1999), [Cair(200^ 
and Cherchneff (2006, 2011, 2012). The post-shock gas consists 
of 1) a narrow layer in the immediate back of the shock front 
where the highly endothermic, collision-induced dissociation of 
molecular hydrogen provides the gas cooling through the reac¬ 
tion H 2 -I- H—+ and 2) a region where the gas cools 
by adiabatic expansion, expands to larger radii, and falls back to 
its initial position owing to the stellar gravity. 

The narrow layer is modelled by using the study of molecular 
shocks in AGB stars by Fox & Wood (1985). We derive the gas 
parameters in the shock front from the Rankine-Hugoniot jump 
conditions for a shocked gas, and assume the shock is weakly 
radiative. The gas cooling is provided by H 2 dissociation, which 
occurs over a space length in the postshock gas given by 


Idiss — "^diss ^ ~ 


1 


^diss ^ **(^ 2 ) 


N,- 


( 1 ) 


where n(H 2 ) is the number density of H 2 in the post-shock gas, 
k;/,jj is the rate of the dissociation reaction of H 2 , V , is the shock 
velocity at the shock formation radius r^, as given in Table T] 
and Njump is the Rankine-Hugoniot velocity shock jump. The 
H 2 dissociation cooling length behind the shock takes values that 
span several orders of magnitude (1 - 10^ cm) depending on the 
shock amplitude (Willacy & Cherchneff 19981. 

Modelling the adiabatic cooling region is based on the study 
by |Berchinger & Chevalier (1985)[ who derived analytical solu¬ 
tions to the hydrodynamic fluid equations that describe the state 
of the gas periodically shocked in the inner wind of a Mira star. 
The inner wind in IK Tau therefore consists of gas layers that 
move forwards and fall back to their initial position over a pul¬ 
sation period, gradually moving outwards owing to shock energy 
dispersion. For the purposes of this study, we assume that the gas 
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Table 2. Gas parameters including the shock velocity Vj, the Mach number M, the gas temperature, and the gas number density in 
the dust formation zone of IK Tau as a function of radius r. Numbers in parenthensis give powers of 10. 


r 


M 

preshock 
no To 

shock front 
n., Tj 

post-shock (/ = 0. 
n, U 

1.0 

32.0 

9.8 

3.6(14) 

2200 

2.1(15) 

42824 

1.2(16) 

6295 

1.5 

26.1 

9.0 

7.8(12) 

1725 

4.4(13) 

28793 

2.2(15) 

4741 

2.0 

22.6 

8.5 

7.2(11) 

1451 

4.1(12) 

21712 

1.9(14) 

3858 

2.5 

20.2 

8.1 

1.4(11) 

1270 

7.7(11) 

17473 

3.4(13) 

3277 

3.0 

18.5 

7.8 

3.9(10) 

1138 

2 .2(11) 

14638 

9.4(12) 

2874 

3.5 

17.1 

7.6 

1.5(10) 

1037 

8 .1(10) 

12630 

3.3(12) 

2568 

4.0 

16.0 

7.4 

6.5(9) 

958 

3.6(10) 

11098 

1.4(12) 

2340 

5.0 

14.3 

7.1 

1.9(9) 

838 

1 .0(10) 

8930 

3.8(11) 

1997 

6.0 

13.1 

6.8 

7.2(8) 

751 

3.9(9) 

7497 

1.3(11) 

1748 

7.0 

12.1 

6.6 

3.4(8) 

684 

1.8(9) 

6460 

6 .0(10) 

1574 

8.0 

11.3 

6.4 

1 .8(8) 

632 

9.9(8) 

5691 

3.1(10) 

1419 

9.0 

10.7 

6.3 

1 .1(8) 

589 

5.9(8) 

5082 

1.8(9) 

1307 

10.0 

10.1 

6.2 

7.1(7) 

553 

3.8(8) 

4597 

1.1(9) 

1209 


Notes. Radius r is in R*, shock velocity in km s temperatures To,s,e in Kelvin, and gas number densities no,j,e in cm 


layers follow strictly periodic motions close to the star and fall 
back to their initial position in the inner wind before the passage 
of the next shock. 

The initial time-averaged conditions of the pre-shock gas at 
each radius follows the formalism developed by Cherchneff et al. 
(1992). For the gas density, the variation as a function of stellar 
radius depends on the extended scale height and is given by 


p(r) = p{r,) X exp 

\ 


G(l-y") 
^o(g)(1 - a ) 



( 2 ) 


where Hoir,) is the static scale height at radius r^, and the param¬ 
eter y is defined as the ratio of the shock velocity and the local 
escape velocity (y = Vjveirs)). The parameter a is the exponent 
of the power law followed by the gas temperature given by 


T(r) = T(r,) X (3) 

where a = -0.6 (IDuari et al. 19991 1. Equation accounts 
for the effect of the shocks on the local gas density 
(ICherchneff et al. 19921 . 

For IK Tau, several shock velocities were tested on the gas 
phase chemistry (20, 25 and 32 km s“'), and the best results 
were achieved for the stronger shock of 32 km s 'at = 1 R*, 
as in the IK Tau study of Duari et al. (1999). Owing to the 
conservation of energy, the shock strength decreases as the 
shock crosses the inner wind according to energy conservation 
and Vs(r) - Vsir^) ^r^jr. The stellar and inner wind param¬ 
eters are summarised in Table We first consider a gas layer 
at ri that is shocked, moves forwards, and falls back to its ini¬ 
tial position owing to the stellar gravity. We run the chemistry 
and grain condensation over two full pulsation periods to check 
for periodicity and then push this gas layer to a larger radius 
r2 - ri + Ar characterised by a lower pre-shock gas density and 
temperature. The outputs at ri, specifically the molecular abun¬ 
dances, as well as the abundance of clusters and dust grains, are 
used as inputs to the calculations at radius r 2 . We run the model 
of the gas phase, the nucleation of dust clusters and their con¬ 
densation in dust grains from = 1 R* to 10 R*. The gas pa¬ 
rameters from 1 R* to 10 R* are listed in Table|^ while the vari¬ 
ation of gas temperature and density over a pulsation phase and 
between 1 R* and 10 R* are illustrated in Figure]^ We direct 
the reader to the papers by Cherchneff et al. (1992), Cherchneff 


(1996), Willacy & Cherchneff (1998), Duari et al. (1999), Cau 
(2002) and Cherchneff (2006, 2011, 2012) for full details of the 
formalism and model we use to describe the inner wind of AGB 
stars. 

This simple analytical model does not represent a fully con¬ 
sistent picture of the shocked upper stellar atmosphere and does 
not aim to reproduce mass-loss rates in AGB stars. We do not at¬ 
tempt to estimate the radiation pressure force acting on newly 
formed dust grains in the gas and derive a mass loss for IK 
Tau. We consider strictly ballistic motions of the gas above the 
photosphere and identify the location of cluster formation and 
dust condensation in the inner wind. Therefore, our model pro¬ 
vides a simplified description compared to more sophisticated 
hydrodynamic models of the inner wind, which take into account 
the dust radiative feedback on the gas (e.g.. [Winters et al. 20001 
IHofner et al. 2003llWoitke 20061) . While such models aim to re¬ 
produce the mass loss rates and optical properties of AGB winds, 
they do not seem to explain the non-equilibrium chemical pro¬ 
cesses that trigger the formation of molecules, dust clusters and 
the growth of dust grains from the gas phase. Indeed, they of¬ 
ten assume over-simplifying mechanisms (e.g., thermodynamic 
equilibrium) for the production of molecules and dust grains. 
Our simplified physical model of the shocked layers above the 
stellar photosphere reproduces to some extent the shocked gas 
conditions, as shown by Bertschinger & Chevalier (1985), and is 
thus sufficient to focus on and understand the chemical processes 
at play in the dust formation zone. 


2.2. Chemical kinetic model and condensation of dust grains 

Photospheric atomic and molecular abundances are calculated 
by assuming TE at 1 R* and by taking the initial elemental com¬ 
position given by the FRUITY database of Cristallo et al. (2011), 
for a 1.5 Mq star on the ZAMS with solar metallicity, that has ex¬ 
perienced three thermal pulses on the AGB. The derived stellar 
mass M* and C/O ratio after these three thermal pulses are 1 Mq 
and 0.6, respectively. The C/O value is close to the assumed C/O 
ratio of 0.75 used in Duari et al. (1999) and Cherchneff (2006). 
The effect of using this lower C/O ratio value as opposed to 0.75 
has been assessed and is minor on all molecular and cluster abun¬ 
dances. We then choose to keep the C/O ratio of 0.75 in order to 
better compare our results with previous models of IK Tau. 
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Table 3. Chemical species and dust clusters included in the IK Tau chemical model.. 


Atoms 




Molecules 




0 

OH 

H 2 O 

02 





C 

CH 

CO 

C 02 

CS 

ocs 

HCN 

CN 

Si 

SiO 

HSiO 

H 2 SiO 

Si 02 

SiS 

Si 2 


N 

NH 

NO 

NH 2 

NH 3 

N 2 



S 

SH 

H 2 S 

SO 

SO 2 

HOSO 



P 

PH 

P 2 

PO 

PN 




Mg 

MgH 

MgO 

MgOH 





Fe 

FeH 

FeO 

Fe 2 





F 

HF 

F 2 






Al 

AlH 

AlO 

AlOH 

AIO 2 

AlCl 

AI 2 

AI 2 O 

Cl 

HCl 

CI 2 

CIO 





Ca 

CaO 

CaOH 

CaCl 





Na 

NaH 

NaCl 






Ti 

TiO 

Ti 02 






K 

KH 

KCl 






H 

H 2 







Dust clusters 

SizOj 

SFOj 

Si303 

31404 

H 2 Si 202 

HSizOj 

H2Si203 

H2Si204 

MgSi 203 

MgSi 204 

MgSi 205 

MgSi 206 H 2 

]VIg 2 Si 204 

Mg 2 Si 205 

Mg 2 Si 206 

Mg2Si207H2 

MgsSFOs 

Mg3Si206 

Mg3Si207 

Mg3Si208H2 

Mg4Si206 

Mg4Si207 

Mg4Si208 


AI 2 O 2 

AI 2 O 3 

A1406 








Phase O 


Fig. 2. Gas temperature and number density variation in the post¬ 
shock gas over one pulsation period as a function of phase for 
radius between 1 and 10 R*. 


2.2.1. Chemistry 


The rich chemistry of the IK Tau inner wind is modelled by a 
network of chemical reactions that includes 100 species. The 
species were chosen according to the following criteria; a) as¬ 
tronomical detection in AGB and supergiant winds, b) associa¬ 
tion with the chemistry of detected molecules, and c) implication 
in the nucleation of dust clusters. The list of species is given in 
Tabled 


All possible types of reactions, effective in a high temper¬ 
ature and density gas, are considered. They include termolecu- 
lar processes that are active at high gas densities, and bimolec- 
ular processes, which include neutral exchange and radiative 
association reactions, and collisional dissociation. The chemi¬ 
cal networks of Duari et al. (1999) and Cherchneff (2006) have 
been revised and extended to include the processes involved in 
the formation of clusters (see below and Sarangi & Cherchneff 
2013, 2015). Reaction rates were taken primarily from the 
National Institute of Standards and Technology Chemical 
Kinetics Database, the 2012 version of the UMIST Database 
for Astrochemistry ( McElroy et al. 20T3]l, specific studies o f ce¬ 
ramic production in flames (e.g., Zachariah & Tsang 1995 i and 
modelling of other circumstellar environments (Cherchneff & 
Dwek 2009, 2010, |Sarangi & Cherchneff 201 3j l. 


Following Cherchneff (2011a, 2012), the implementation of 
rates for reverse reactions was changed to accommodate avail¬ 
able new rates measured in combustion and aerosol chemistry. 
The rate of a specific reverse process was not calculated from 
the equilibrium constant by assuming detailed balance as in 
Duari et al. (1999) and Cherchneff (2006), but directly en¬ 
tered in the chemical network with its available measured, cal¬ 
culated, or estimated value. This new, purely kinetic approach 
has been tested for the case of the carbon star IRC-110216 
and provides better results when compared to observational 
data. Specifically, the formation of water and other relevant 
molecules in carbon stars was successfully explained by this re¬ 
vised shock-induced chemistry applied to the inner stellar wind 
(ICherchneff 201 lallCherchneff 20121 . 

No ionic species and related chemical reactions are consid¬ 
ered because we assume there are no ionisation processes of 
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Table 4. Chemical reactions that initiate the formation 
of silaformyl and the silanone dimerisation, adapted from 
Zachariah & Tsang (1995), and the prevalent reactions leading 
to the formation of enstatite dimer, adapted from Goumans & 
Bromley (2012). 


R1 

SiO 

+ 

H 


HSiO 



R2 

SiO 

+ 

H 2 


HSiO 

+ 

H 

R3 

HSiO 

+ 

H 


H2SiO 



R4 

HSiO 

+ 

H 2 


H2SiO 


H 

R5 

HjSiO 

+ 

H2SiO 


H2Si202 

+ 

H 2 

R6 

HjSLOz 

+ 

H 2 O 


H2Si203 


H 2 

R7 

HjSLOj 

+ 

Mg 


MgSLOs 

+ 

H 2 

R8 

MgSLOs 

-f- 

H 2 O 


MgSi204 

+ 

H 2 

R9 

MgSi 204 


Mg 


Mg 2 Si 204 



RIO 

Mg2Si204 

+ 

H 2 O 


Mg2Si205 

+ 

H 2 

Rll 

Mg 2 Si 205 

+ 

H 2 O 


Mg2Si206 

+ 

H 2 


atomic and molecular species induced by the ultraviolet stellar 
radiation field, which is low for the effective temperature of IK 
Tau. Furthermore, we assume there is no penetration of the in¬ 
terstellar radiation field or cosmic rays in these regions close to 
the stellar surface. 

2.2.2. Dust clusters 

Together with the formation of molecules, we model the synthe¬ 
sis of molecular clusters that are precursors to dust grains. We 
consider two types of prevalent condensates in the inner wind, 
silicates and alumina, and consider as dust clusters the dimers 
of forsterite (Mg 2 Si 04 ) and enstatite (MgSiOs), and the dimer 
of AI 2 O 3 . The formation routes of silicate dimers of forsterite 
and enstatite stoichometry were studied by Goumans & Bromley 
(2012). By identifying the structures of particularly energetically 
stable small clusters, their thermodynamic properties were esti¬ 
mated and used to highlight potential cluster formation routes at 
equilibrium. For the gas parameters illustrated in Figure the 
exergonocity of reactions involving clusters is assessed by us¬ 
ing the thermodynamic data estimated by Goumans & Bromley 
(2012) for the clusters listed in Table and reaction rates for 
both formation and destruction processes are estimated on the 
basis of known rates for similar reaction types and the likelihood 
of the reactions to proceed. The synthesis of silicate clusters is 
described by the recurrent addition of oxygen to the cluster lat¬ 
tice through reaction with H 2 O, and the inclusion of one Mg 
atom after the oxidation step. 

This suite of reactions leads to the formation of enstatite 
(MgSiOs) and forsterite (Mg 2 Si 04 ) dimers. In the study by 
Goumans & Bromley, the growth to silicate small clusters fol¬ 
lows a first important step, the dimerisation of SiO, which acts 
as a bottleneck to the growth mechanism. However, the dimeri¬ 
sation of SiO is a slow process. No measured rate value ex¬ 
ists, but the theoretical study by Zachariah & Tsang (1993) pro¬ 
vides estimates of SiO addition rates up to the SiO pentamer, 
and their reverse dissociative paths, and these rates are used in 
this study. Because of the overwhelming abundance of H and 
H 2 in the inner wind, another pathway to the formation of small 
Mg/Si/O clusters is included that involves the dimerisation of the 
silaformyl radical, HSiO, followed by the formation of silanone, 
H 2 SiO, and its dimerisation. These chemical routes were studied 
by Zachariah & Tsang (1995) in the context of the nucleation of 
silica in a silane-rich flame. We list in Table|4]the chemical reac¬ 
tions involved in the nucleation of the MgSi 203 molecular clus¬ 
ter (R1-R7). The later growth of small silicate clusters resumes 


with the mechanism derived by Goumans & Bromley, and the 
oxygen addition to MgSi 203 through reaction with the prevalent 
O-bearing species, H 2 O, according to Reactions R8-R11. 

The formation of alumina dimers, (Al 203 ) 2 , is described 
as follows: the formation of AI 2 O 3 results from the trimolec- 
ular dimerisation of AlO, followed by the oxidation of (A 10)2 
through reaction with O-bearing species, essentially H 2 O. This 
choice was based on similarities in the structure of Si 203 
and AI 2 O 3 as both molecules have kite-shaped, global mini¬ 
mum structures (IPatzer et al. 19991 Li & Cheng 2012 1 . Because 
no measurements appear to exist on the gas-phase formation 
routes to AI 2 O 3 , we have assumed chemical routes and related 
rates similar to those for the formation of Si 203 and H 2 Si 203 . 


2.2.3. Condensation of dust grains 

Simultaneously to the formation of molecules and dust clusters 
from the gas phase by non-equilibrium chemistry, we model 
the condensation of dust grains following Plane (2013) and 
Sarangi & Cherchneff (2015). We choose the forsterite and alu¬ 
mina dimers, (Mg 2 Si 04)2 and (Al 203 ) 2 , respectively, as sta¬ 
ble condensation nuclei, and consider both the coalescence and 
coagulation of these dimers to form dust grains. The formal¬ 
ism is based on a Brownian formalism, which accounts for the 
collision, scattering, coalescence and coagulation of the grains 
through Brownian motion. We use the equation formulated by 
Jacobson (2005), whereby the variation with time of the num¬ 
ber density of a grain of specific volume v is described by the 
integro-differential coagulation equation given by 


drivit) 

dt 


V 


Pv-v^vnv-v'nv'dv' 


00 


dv'. 


(4) 


Here, t is the time, v' and (v-v') are the volumes of the two coag¬ 
ulating particles, n^ is the number density of grains with volume 
V, and /3vy is the rate coefficient of coagulation between particles 
with volume v and v'. The shocked gas is characterised by a free- 
molecular regime, defined by the grain mean-free path in the gas 
being much larger than the grain size, where Brownian diffusion 
prevails in the coagulation process and controls the rate of coag¬ 
ulation for particles i and j. Equationj^was solved by using 
a semi-implicit volume conserved model, where grains are as¬ 
sumed spherical (ISaunders & Plane 2006t . For both alumina and 
silicate, grains are assigned to discrete bins, following a volume 
ratio distribution given by v„ = Y'^^vq, where v„ is the volume 
of the bin, y is a constant defined as the ratio of the volumes 
of adjacent bins, and vq is the volume of the first bin, determined 
by the size of the gas-phase precursor. At each time step, the vol¬ 
ume of the first bin, vq corresponds to the radius qq of the stable, 
largest clusters produced from chemical kinetics in the nucle¬ 


ation phase (for alumina and silicate, uq = 3.45 A and 3.3 A, re¬ 
spectively - see Sarangi & Cherchneff 2015) . The grain size dis¬ 
tributions are then derived for silicates of forsterite and enstatite 
stoichiometry, and alumina. We direct the reader to Plane (2013) 
(and references therein) and Sarangi & Cherchneff (2015) for a 
detailed description of the condensation formalism, the expres¬ 
sions for the rates of Brownian diffusion and coalescence, and 
their variation with gas parameters. 

Our formalism for grain growth does not consider deposition 
of atoms and molecules on the surface of small grains. Once 
the radius of grains reaches several tens of A, the growth of 
the particles proceeds through surface deposition, which is a 
fast process, and coagulation and coalescence, as observed in 
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flames jRoth 2007t . For deposition, the sticking coefficient S of 
an atom or molecule on the surface of grains depends on the 
chemical composition of the accreting species and the grain sur¬ 
face, and the surface and gas temperature. The quantity may vary 
over orders of magnitude depending on the grain material, with 
value range of ~ 10“"^ - 1 (e.g., |Leitch-Devlin & Williams 1985] 
IBrune et al. 19931 [Nagahara & Ozawa 1996| l. Existing studies 
on dust formation in circumstellar environments often assume a 
growth process through surface deposition only, and an S value 
equal to unity or in the range 0.1-1 (e.g., Gail & Sedlmayr 1999| 
IBladh et al. 20151) . In view of these uncertainties, we prefer to 
not consider surface deposition in the formalism at this stage of 
our study, and discuss the potential effects of accretion on the 


present results in § 3.2 


3. Results 

The derivation of molecular and cluster abundances is as fol¬ 
lows: we assume that the initial composition of the pre-shocked 
photosphere is given by TE calculations for the parameters in 
Table [T] We derive the abundances of molecules and clusters in 
the immediate “fast chemistry” post-shock gas and the“adiabatic 
trajectory” of a gas parcel (see Willacy & Cherchneff 199^ by 
solving a system of 100 master equations described by a chemi¬ 
cal network that simultaneously includes the formation of chem¬ 
ical species, and alumina and silicate clusters. When the clusters 
have reached the size of dimers of alumina and forsterite, they 
start to grow by coalescence and coagulation. As explained in 
the parcel of gas is then lifted from radius r to radius r', and 
the calculated abundances for the shocked gas at ri are then used 
as initial pre-shock conditions for the gas parcel at r 2 - + Ar, 

with Ar - 0.5 R*. The radius ri was then varied from 1 R* to 
10 R*. We first discuss the results of the gas phase in § |3.1| and 
then present results on dust clusters and dust grains in §|3.2| 




3.1. Gas phase 

Molecular abundances versus pulsation phase at two radii in 
the inner wind are shown in Figure At 1 R*, equivalent to 
the shocked stellar photosphere position, the shock destroys all 
molecules in the H 2 dissociation zone, and species gradually 
reform in the cooling region until p ^ 0.2 - 0.4, depending 
on molecules. Beyond this phase, some molecules are fully re¬ 
formed with almost constant abundances, such as CO, SiS, and 
PN, while others are still chemically active, as for example, SO, 
and AlO. At larger radius, the early destruction of molecules in 
the H 2 dissociation zone still occurs, but it is less efficient. Some 
molecules maintain almost similar abundances to those at 1 R*, 
e.g., SiS, SiO, and PN, while others are still chemically active 
(e.g., AlO). The prevalent molecules formed in the inner wind 
are illustrated in Figure]^ where abundances with respect to to¬ 
tal gas as a function of radius are shown, at a phase (p - 1. As 
previously mentioned, some molecules form in the post-shock 
gas of the shocked photosphere and maintain abundances al¬ 
most constant in the inner wind. These species include CO, PN, 
HCl, TiO, and SiS. The rather constant abundances reflect the 
fact that those species do not participate in complex chemistries 
and in the formation of dust clusters in the post-shock gas. Once 
they form in the shocked gas at 1 R*, they stay unchanged until 
the chemistry is frozen out at the wind acceleration radius after 
r > 6 R* (see § |3.2.2| ). 

Conversely, other molecules actively partake in the chem¬ 
istry of the gas phase, including that of dust cluster formation. 
This is the case for H 2 S, SO, AlOH, SiO, and the carbon-bearing 


Fig. 3. Molecular abundances with respect to total gas number 
density n^a , as a function of pulsation phase for two specific radii 
in the IK Tau inner wind. Top: 1 R*; Bottom: 3 R*. 



Fig. 4. Molecular abundances with respect to total gas number 
density n^aj as a function of radius in the inner wind of IK Tau. 
Values correspond to = 1 in the gas layer and the grey shaded 
area corresponds to the region where the wind is already accel¬ 
erated - see § |3.2.2| 
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CO SiO PN SO NaCI PO 


Fig. 5. Modelled and observed molecular abundances in the in¬ 
ner wind of IK Tau. The modelled values are for (p - and 
r = 6 R* and 8 R*. The bars for observed data indicate the 
value range derived from the different sets of observations listed 
in Tabled 


S H- CN ^ CS -I- H, 


(11) 


and 

CO-H SH ^ OCS-H H, (12) 

followed by 

OCS-H H ^ CS -I- OH. (13) 


Therefore, the CS, CN, SH and H 2 S molecules are chemically 
linked. Inspection of Figurej^indicates a dip in the CS and HCN 
abundances between 2 R* and 4 R*. In that radius range, the 
formation of silicate dust clusters that is initiated by the forma¬ 
tion of HSiO and H 2 SiO depletes molecular hydrogen to release 
atomic hydrogen in the process. This excess favours the reverse 
of reaction]^ producing copious amounts of sulphur atoms and 
decreasing the SH and H 2 S reservoirs. Therefore, the CS abun¬ 
dance drops as reactions 12 and 13 become less efficient. The 
CN radical initially forms in the H 2 dissociation region after the 
shock front from the reactions 


N-hCO^CNh-O, (14) 

and 


species CO 2 , HCN and CS. While SiO and AlOH, through AlO, 
are directly associated with the nucleation of silicates and alu¬ 
mina, respectively, other molecules are linked by the shock- 
sensitive chemistry, like H 2 S and SO. The formation of H 2 S oc¬ 
curs close to the star through the reaction 

SH -H H 2 ^ H 2 S -H H, (5) 

where SH forms in the hot post-shock gas through the reaction 

S -H H 2 ^ SH H- H, (6) 

which has a high activation energy. At larger radii and for weaker 
shocks and lower post-shock temperatures, the formation of SO 
proceeds essentially through 

S -H OH ^ SO -H H, (7) 


N -H CS ^ CN H- S. 


(15) 


The formation of HCN occurs in all post-shock regions of 
the inner wind following the reaction 


CN -H H 2 ^ HCN -H H 


(16) 


Both CN and HCN abundances are regulated by reaction 16 
and its reverse process, and the formation of HCN is then con¬ 
trolled by the H 2 molecule in the inner wind. Therefore, the drop 
in the HCN abundance ensues from the depletion of H 2 due to 
the nucleation of silicates between 2 R* and 4 R*. 

Finally, the molecule CO 2 was detected with the 
Infrared Space Observatory, ISO, in two Mira stars, T Cep 
(lYamamura et al. 1999t and R Cas (IMarkwick et al. 2000t . 
but has not been detected in IK Tau to date. We predict the 
formation of CO 2 by the reaction 


for which the only rate available is at 300K, and 

SH -I O ^ SO -I- H. (8) 

The rates for reactions |7] andhave no activation energy barrier 
and can easily proceed in the inner wind to form SO. However, 
the reverse reactions have a high activation energy (~ 11000 K 
and 20000 K for the reverse of reactions [7] and respectively), 
and thus proceed close to the star and inhibit the formation of 
SO at small radii. Through reaction SO formation depletes 
the SH reservoir and then hampers the formation of H 2 S, whose 
abundance starts dropping around 3 R*. The production of SO 
triggers the formation of SO 2 from the gas phase following the 
reaction 

SO -H OH ^ SO 2 H- H, (9) 

which is barrierless and can proceed as soon as the SO abun¬ 
dance builds up. However, we see from Figure that the SO 2 
abundance remains always low in the inner wind, and much 
lower than that of SO by almost two orders of magnitude. 

The C-bearing species, HCN and CS, form from a complex 
chemistry that is linked to the SH/H 2 S chemistry. Basically, the 
dominant formation processes for CS are 

Ch-SH^CS-hH, (10) 


CO-H OH ^ CO 2 -H H, (17) 

which is active in O-rich AGB winds dCherchneff 20061 1. This 
reaction goes backwards when the synthesis of silicate clusters 
occurs and replenishes the gas in atomic H at the expense of H 2 
in the 3 - 6 R* region. A slight drop in CO 2 abundance then 
occurs, as shown in Figure]^ 

The modelled molecular abundances for IK Tau are sum¬ 
marised in Table 1^ where abundances in the stellar photosphere 
at TE are listed, and the abundances in the shocked wind are 
given at 1 R* and 6 R*. These values are compared to abun¬ 
dances derived from the modelling of observational data of the 
IK Tau inner wind. We see that two distinct groups of molecules 
exist. The hrst group includes species that are present at TE in 
the photosphere, are destroyed by shocks, but reform with al¬ 
most their TE abundance values (e.g., H 2 O, SiO, SiS, SO, CO 2 , 
TiO, and HCl). The second group includes species that are not 
present in the photosphere and only form in the inner wind as a 
result of shock chemistry (e.g., HCN, CS, PN, AlOH, NaCl). As 
for the agreement between modelled and observed abundance 
values, we consider that both sets of values agree well when the 
discrepancy is within a factor of 10. In Eigure we plot the 
modelled abundances at 6 R* and 8 R* and the values derived 
from observations. We see that both modelled and observed sets 
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Table 5. Predicted abundances (with respect to H 2 ) of selected molecules. TE abundances are given at r^ and modelled abundances 
are listed for r^ and 6 R*. Most recent abundances (with respect to H 2 ) of molecules observed in the inner wind of IK Tau and other 
objects are listed, with the corresponding references. The double horizontal line separates data available for IK tau (upper part) from 
data on other objects (lower part). 


Species 

TE at r. 

Modelled at r^ 

Modelled at 6 R* 

Observed 

Reference 

CO 

1.0 x 10^’ 

7.3 X lO^** 

7.4 X 10-“ 

2.0 X 10-'' 

Decin et al. 2010bl 





3.0 X 10-'' 

Kim etal. 20101 

H 2 O 

2.1 X lO-** 

1.1 X 10-* 

1.8 X 10-* 

6.6 X 10-5 

Decin et al. 2010bl 

SiO 

9.5 X 10-" 

6.4 X 10^" 

1.9 X 10-5 

8.0 X 10-“ 

Decin et al. 2010bl 


- 

- 

- 

5.1 X 10-“ - 1.3 X 10-5 

Kim etal. 20101 

SiS 

4.0 X 10-** 

2.2 X 10-7 

2.7 X 10-7 

1.1 X 10-5 

Decin et al. 2010al 


- 

- 

- 

3.1 X 10-’ - 1.3 X 10-“ 

Kim etal. 20101 


- 

- 

- 

1.0 X 10-’ 

Schoier et al. 20071 

SO 

4.3 X 10-** 

2.1 X 10-‘^ 

7.3 X 10-** 

4.0 X 10-' 

Decin et al. 2010al 



- 

- 

2.7 X 10-’ - 7.8 X 10-’ 

Kim etal. 20101 

S 02 

1.3 X 10^" 

1.8 X 10-“* 

4.3 X 10-'^ 

2.0 X 10-“ 

Decin et al. 2010aJ 


- 

- 

- 

4.2 X 10-“ - 1.4x 10-5 

Kim etal. 20101 


- 

- 

- 

1.7 X 10-**-3.7 X 10-’“ 

Yamamura et al. 19991 

HCN 

6.3 X 10^'^ 

7.7 X 10-' 

5.0 X 10-** 

1.3 X 10-' 

Schoier et al. 201.3j 


- 

- 

- 

4.3 X 10-’ - 1.4x 10-“ 

Kim etal. 20101 


- 

- 

- 

4.4 X 10-’ 

Decin et al. 2010bl 

CS 

1.9 X 10-" 

2.5 X 10-'’ 

5.7 X 10-' 

8.0 X lO-** 

Decin et al. 2010bJ 


- 

- 

- 

8.1 X 10-**-3.0 X 10-’ 

Kim etal. 20101 

PN 

4.0 X 10-'“ 

5.7 X 10-’ 

5.8 X 10-'' 

3.0 X 10-’ 

De Beck et al. 20131 

PO 

1.0 X 10-' 

3.3 X 10-"’ 

2.0 X 10-"’ 

5.4 X 10-’'-6.0x 10-’ 

De Beck et al. 20131 

NaCl 

3.7 X 10^'^ 

7.8 X 10-"’ 

1.0 X 10-** 

4.0 X 10-** 

Milam et al. 20071 

NH 3 

1.1 X 10-" 

2.8 X 10-« 

1.5 X 10-*** 

1 - 3 X 10-“ 

Menten et al. 20101 

n 

p 

4.1 X 10-** 

2.5 X 10-** 

4.4 X 10-** 

5.9 X 10-** * 

Markwick et al. 20001 


- 

- 

- 

2 X 10-“ - 4 X 10-** “ 

Yamamura et al. 19991 

AlO 

8.0 X lO^'’ 

2.6 X 10-* 

1.6 X 10-'“ 

3.4 X 10-**“ 

Kaminski et al. 2013al 

AlOH 

4.5 X 10^® 

4.5 X 10-'’ 

1.5 X 10-“ 

1.0 X 10-’ “ 

Tenenbaum & Ziurys 2010| 

TiO”^ 

2.4 X 10-’ 

1.8 X 10-’ 

1.8 X 10-’ 

5.0 X 10-“ “ 

Kaminski et al. 2013bl 

Ti 02 

1.3 X 10-‘2 

6.2 X 10-’“ 

4.8 X lO-"* 

6.3 X 10-'“ “ 

Kaminski et al. 201.3bl 


- 

- 

- 

3.8 X 10-** “ 

De Beck et al. 20151 

HCl 

3.5 X 10-’ 

3.9 X 10-’ 

3.7 X 10-’ 

1.5 X 10-*“ 

Yamamura et al. 20001 

AlCl 

1.5 X 10-'“ 

3.8 X 10-‘2 

2.2 X 10-'“ 

< 5 X 10-* “ 

Tenenbaum & Ziurys 2010| 

SH 

2.8 X 10-^ 

1.3 X 10-’ 

2.3 X 10-** 

2.0 X 10-’ “ 

Yamamura et al. 20001 

HjS 

8.1 X 10-® 

2.8 X 10-5 

1.0 X 10-** 

< 2 X 10-“ “ 

Yamamura et al. 20001 


Notes. SO 2 & CO 2 abundances derived from data of T Cephei ® The [C02]/[H20] ratio as measured in R Cas was used to estimate CO 2 
abundance in IK Tau Measured in the supergiant VY Canis Majoris Also detected in the inner wind of the S star NP Aurigae by Smolders et 
al. (2012) Measured at 1.1 R* in the S starR Andromedae. 


of data agree well for 9 species out of 11 molecules detected. On 
the other hand, the modelled abundance values for SO 2 and PO 
fail to agree with observations. 

For SO 2 , we find a low abundance of ~ 2 x 10“® at 6 R*, 
which indicates that the molecule does not efficiently form in 
the gas phase from the reaction of SO with OH. Yamamura et al. 
(1999) detected the 7.3 /tm band of SO 2 in the ISO/SWS spec¬ 
tra of the Mira star T Cep, and derived an SO 2 abundance in the 
range 2 x 10“* - 4 x 10“^ in a shell located at 4 - 5 R*. Most of 
the IR band emission comes from the inner layers, and we notice 
that our modelled SO 2 abundance at 4 R* is ~ 4 x 10“^, a value 
consistent with the T Cep observations. Decin et al. (2010a) de¬ 
rive a high abundance of 1 x 10 ® from fitting emission lines 
originating in the wind between 50 R* and 200 R*. This region 
is located outside the inner wind and overlaps with the outer 
envelope, defined as the region where the interstellar ultravi¬ 
olet radiation field starts dissociating molecules and triggers a 
cold, ion-rich chemistry. Combined with the penetration of cos¬ 
mic rays, the region starts around 100 R*. The high abundance 
of SO 2 derived from observations may not reflect a high molec¬ 


ular abundance in the inner wind at r < 10 R*, but a possible 
high molecular abundance at the onset of the outer envelope. 

The region between 10 R* and 100 R*, labelled “the in¬ 
termediate envelope” in Figure [T] is characterised by cool gas 
temperatures (~ 100 K), and the presence of parent molecules 
and dust formed in the inner wind. It has long been hypothe¬ 
sised that some molecules commonly detected in AGB winds, 
e.g., methane, CH 4 in carbon stars and ammonia, NH 3 , in 
AGBs, could originate from this region (lOmont 19881 1. Indeed, 
for IK Tau, we form very small amounts of NH 3 from the 
gas phase in the inner wind, as seen from Table while a 
large abundance has been derived from recent Herschel obser¬ 
vations (IMenten et al. 201 Oi l. These data indicate a formation re¬ 
gion with gas number densities of ~ 10 ® cm~^ and gas tempera¬ 
tures in the range lO-lOOK, typical of the intermediate envelope. 
Therefore, molecules that cannot form in the shocked inner wind 
from the gas phase could well be synthesised on the surface of 
dust grains in the intermediate envelope. 

A different scenario applies to PO. The phosphorous chem¬ 
istry is very poorly characterised and we have simply used the 
isovalence between nitrogen and phosphorous to derive some 
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Fig. 6. Modelled abundances of Al-bearing species and AI 2 O 3 
clusters at r = 1 R*. 


chemical rates for P-bearing species. While the modelled abun¬ 
dance for PN slightly exceeds but agrees well with the observed 
value, the modelled PO abundance is ~ 3 orders of magnitude 
smaller than the value derived from observations by de Beck et 
al. (2013). In the phosphorous chemistry, the dominant forma¬ 
tion process for PN is the reaction 

POh-N^PNh-O, (18) 

while PO essentially forms from the reaction 
P -H OH ^ PO -I- H, (19) 

and 

PN -H OH ^ PO -I- NH. (20) 

The reaction rates of reactions [TSl-ISOl have been estimated from 
the N/P isovalence and have not been measured or theoretically 
assessed. However, the formation of PO is linked to that of PN. 
Then, the slight overestimate of the depletion of P atoms in PN 
as shown by our modelled value may result in lowering the PO 
abundance. We thus attribute the discrepancy between modelled 
and observed PO abundances to a poorly characterised chemistry 
but note that the isovalence principal gives a satisfactory match 
for PN. 

3.2. Dust 



a(A) 

Fig. 7. Grain size distributions for alumina dust formed in the 
shocked gas layer. Top; as a function of phase in one oscillation 
over one pulsation period at 1 R*; Bottom; the shocked layer has 
moved from 1 R* to 2 R*. Distributions are shown for different 
pulsation phases (p (with p - 1 - 6 ) at radius r = 1 - 2 R*. 


at IR* may not survive the passage of the next shock after one 
pulsation period. Any dust formed at IR* may thus not accu¬ 
mulate and grow over several pulsation periods in the shocked 
photosphere. The situation is different at larger radii, where the 
post-shock gas temperature drastically decreases. 


The formation of the solid phase is coupled to the gas-phase 
chemistry and calculated at each integration time step in the 
post-shock gas. We present results on grain size distributions for 
both alumina and silicate grains, and list the derived dust-to-gas 
mass ratio in Table The build-up of dust grains occurs over a 
number of pulsations that are necessary for the gas layer to move 
from r R* to (r - 1 - 6r) R*. We calculate this number of pulsations 
by assuming a local gas drift velocity. The drift velocity is small 
close to the star owing to stellar gravity, and increases further 
away from the stellar surface. We choose drift velocity values 
derived from pulsation models, and assume a drift velocity of 
0.5 km s ' for r < 2 R* and of 1.5 km s ' for r > 3 R*. 

Owing to the high gas temperatures characterising the post¬ 
shock gas at 1 R* (see Figure]^, it is reasonable to consider that 
the dust clusters and small grains that form in one oscillation 


3.2.1. Alumina 


The alumina clusters chiefly form in the post-shock gas at 1 R*, 
and their abundances in the first oscillation are illustrated in 
Figure The formation of AI 2 O 3 involves the dimerisation of 
AlO through termolecular recombination while AlO forms from 
the reaction of Al with OH in the first place. The dimerisation 
of AI 2 O 3 is also described by termolecular recombination to 
form the AI 2 O 3 dimer, (Al 203 ) 2 , which has a binding energy 
of ~ 15.7 eV (Li & Cheng 2012 1 , and can thus withstand the 


high temperatures of the post-shock gas. This dimer is consid¬ 
ered as the seed that will further condense into alumina grains. 
Over one oscillation resulting from the first shock at 1 R*, only a 
small quantity of alumina dimers forms, reaching an abundance 
of 1 X lO^'" at(f> - 1. A large fraction of Al stays in the form of 


10 















Gobrecht et al.: Dust formation in IK Tau 


molecular rings (e.g., AI 2 O 2 and AI 2 O 3 ). Because the dimerisa- 
tion of AlO and AI 2 O 3 is termolecular and thus dependent on the 
background pressure, cluster formation occurs over a reduced 
range of radii in the inner wind and the formation of AI 2 O 3 clus¬ 
ter ceases after 2 R*. 

The grain size distributions for alumina in one oscillation at 
1 R* is presented in Figure [7] as a function of pulsation phase. 
Because of the high post-shock gas temperatures after the shock 
passage and the extremely low abundance of alumina dimers, no 
grains form at 0 < 0 < 0.5. Very small grains start forming after 
4> - 0.5, but coagulation starts to be efficient for dust growth in 
the phase range (j) - 0.8 - 1 , for which gas temperatures span 
the range ~ 1600 - 2000 K and the alumina cluster abundance 
reaches a maximum. The distribution peaks at grains in the size 
range 70 - 200 A, despite the small fraction of AI 2 O 3 clusters 
available. However, the post-shock gas densities are very high 
at 1 R*, as seen from Table [T] and the clusters coagulate and 
condense with high efficiency and lead to the formation of AI 2 O 3 
grains of small size over one oscillation and pulsation period. 
Under the assumption that small alumina grains may not survive 
the passage of periodic shocks, the results show that there exists 
a steady population of small size {a ~ 0.02 /rm) alumina grains 
located in the shocked photosphere at 1 R*. If we consider that 
grain growth may also proceed through surface deposition when 
the grain radius has reached several tens of A (at 0 = 0.8 and 
after), the grain size will increase, but the growth will be time- 
limited and proceed only between ~ 0.8 < 0 < 1 . 

We now consider the extreme case where alumina grains sur¬ 
vive the periodic shocks at each shock passage, while molecules 
are destroyed and reform. We may then investigate grain growth 
over several pulsation periods. For a gas drift velocity of 0.5 km 
s ', the gas layer at 1 R* takes 12 pulsations to reach 2 R*. We 
thus consider the build-up of dust grains over a number of 12 
pulsations for the standard gas conditions given in Table and 
show the resulting alumina grain size distributions in Figure 
The grain size distributions are skewed towards large grain sizes. 
At 2 R*, the distribution comprises two dust grain components 
peaking at ~ 300 A and 0.3 //m. All small grains produced in the 
first oscillation at the shocked photosphere have effectively con¬ 
densed in larger grains after 12 pulsations. The distribution does 
not change at r > 2 R*, owing to a shortage of growing agents as 
all alumina clusters have been included in dust grains. The mass 
of alumina dust produced remains modest and the dust-to-gas 
mass ratio varies from 3.3 x 10“^ at 1 R* over one oscillation to 
2.8 X 10“^ at 2 R* when grain survival is assumed, as seen from 
Table | 6 ] 


3.2.2. Silicates 

Small clusters of silicates form further out of the shocked, stel¬ 
lar photosphere, as seen in Figure where abundances of key 
molecules involved in the nucleation of small silicate clusters 
are shown as a function of radius. The silaformyl radical, HSiO, 
and silanone, H 2 SiO, start to efficiently form at r > 1.5 R* be¬ 
cause the reverse processes of reactions R1-R4 in Table are 
favoured and the destruction of these nucleating species ensues 
at shorter radii. The synthesis of dimers of forsterite stoichiom¬ 
etry starts at r > 3 R*, and the abundances of the prevalent clus¬ 
ters involved in the nucleation phase are shown in Figure]^ for 
the gas at r = 3.5 R*. The gas temperature at which all inter¬ 
mediate clusters are effectively converted into silicate dimers is 
~ 1200 K in the post-shock gas, where the gas number density 
is ~ 4.5 X lO'" cm“^. 




Fig. 8. Top; Abundances with respect to total gas number den¬ 
sity ngfli of the silicate clusters and key molecules as a function 
of radius; Bottom; Abundances with respect to total gas num¬ 
ber density n^^j of the silicate clusters as a function of pulsation 
phase at r = 3.5 R*. 


The build up of silicate dust takes place from 3.5 R* to 
larger radii, over several stellar pulsations, as we assume that the 
small silicate grains formed in one oscillation at 3.5 R* with¬ 
stand the next shock passage. As mentioned in § 3.2 for a gas 
drift velocity of 1.5 km s“', it takes 2 pulsations to travel a dis¬ 
tance of 0.5 R*. Therefore, 10 pulsations are necessary to move 
from 3.5 R* to 6 R*, where the growth of silicate grains slows 
down and almost reach completeness at 8 R*. Figure shows 
the silicate grain size distributions at various radii and pulsa¬ 
tion numbers. We see that the size distribution has a tail of large 
grains that peaks around 50 A at 3.5 R* and reaches ~ 200 A at 
r > 6 R*. The dust-to-gas mass ratio listed in Table reaches 
1.0 X 10“^ at 6 R*, and the ratio keeps increasing out to 10 R*, 
where it reaches 2.3 x 10^^. These values agree with dust-to- 
gas mass ratios derived from observations of O-rich Mira stars. 
Based on the derived gas-to-dust mass ratios, we estimate that 
~ 22 % of the elemental silicon abundance is depleted in silicate 
grains at 6 R*. 

Because the growth of silicate grains takes place between 
3 R* and 6 R* and slows down starting at 6 R*, the wind accel¬ 
eration, which is triggered by radiation pressure on dust grains, 
should start developing as soon as enough dust grains form in the 
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Table 6. Derived dust-to-gas mass ratio for alumina and silicate dust grains as a function of radius for the gas number densities 
given in Table and the gas number densities increased by a factor of 10 (labelled“xlO”)- The radius range where strong wind 
acceleration starts is indicated in italic. 


r(R*) 


Alumina 



Silicate 


1 Osc. 

Accumulated 

1 Osc. (xlO) 

Accumulated (xlO) 

Accumulated 

Accumulated (xlO) 

1.0 

3.3 X 10-'’ 

2.1 X 10^" 

6.7 X 10^" 

4.3 X 10-“* 

- 

- 

1.5 

2.6 X 10^''' 

2.7 X 10^^ 

5.4 X 10^'’ 

5.3 X 10-^ 

- 

- 

2.0 

6.7 X 

3.3 X 10-^ 

3.6 X 10-*’ 

6.5 X 10-* 

- 

- 

3.5 



” 


1.3 X 10-* 

4.7 X lO^* 

4.0 





4.6 X lO-'* 

1.1 X 10-3 

5.0 



” 


7.3 X lO-'* 

3.0 X 10-3 

6.0 



” 


1.0x10-^ 

5.1 X 10-^ 

7.0 





1.4x10-^ 

6.1 X 10-^ 

8.0 



” 


1.7x10-^ 

6.7x10-^ 

9.0 



” 

” 

2.0x10-^ 

7.5 X 10-^ 

10.0 




” 

2.3 X 10-^ 

8.0 X 10-^ 


Notes. “1 Osc.” corresponds to the case where the dust forms over one oscillation but is destroyed in the next shock.“Accumulated” corresponds 
to the case where the dust forms over one oscillation, is not destroyed by the next shock and accumulates over several pulsation periods - see text. 



Fig. 9. Grain size distributions for silicate dust formed in a 
shocked gas layer that has moved from 3.5 R* to 10 R*. 


radius range 3.5 - 4 R*, and increase to reach full regime when 
the grain sizes are large enough. As for alumina, the deposition 
of chemical species on the grain surface will also proceed along 
with coalescence and coagulation as soon as the grain radius ex¬ 
ceeds several tens of Angstroms. This may happen as near as 
3.5 - 4 R* as shown in Figure]^ where some dust grains reach 


a radius of ~ 50 A at 3.5 R* and ~ 100 A at 4 R*. We expect 
that deposition will push the final grain size distribution to radii 
a > 200 A, and increase the final dust-to-gas mass ratio, but 
cannot yet quantify the increase. From Figure]^ we note that the 


modelled SiO abundance derived at 6 R* 


twice as large as 


the value derived from observations. The modelled SiO value is 
then consistent with the potential growth of dust grains through 
gas deposition along with ongoing coalescence and coagulation. 
The prevalent wind acceleration zone may then encompass the 
radius range 3.5 -6 R*. This acceleration will result in the freez¬ 
ing of molecular abundances after r > 6 R*. This is the underly¬ 
ing reason for comparing modelled abundances at 6 - 8 R* with 
values derived from astronomical data in § 3.1 and Figure]^ 


3.2.3. A high gas number density case 

Despite the fact that the present standard model for IK Tau pro¬ 
vides satisfactory dust-to-gas mass ratios for silicates compared 
to values derived from observations, the grains in which most of 
the silicate mass is concentrated are of medium size (i.e., ~ 200 
A). We thus test the impact of the gas number density on the 
synthesis of clusters and dust grains by increasing the gas num¬ 
ber density at the photosphere by a factor of ten and keeping all 
photospheric parameters and the shock characterisation as be¬ 
fore. This increment in number density is certainly too high to 
properly describe the IK Tau wind at the base of the interme¬ 
diate envelope, defined as where the wind is fully accelerated. 
Assuming the mass loss rate and wind terminal velocity given 
in Table [T] we derive from the definition of mass loss rate and 
conservation of mass. 


n(r) = 


M 

4nmHfigasVexpr^ 


1.09 X 10'° 

yfl 


( 21 ) 


where mn is the mass of the hydrogen atom, Ugas the gas mean 
molecular weight at the acceleration radius (taken as 2 a.m.u. be¬ 
cause the gas is essentially in the form of H 2 ), Vexp the terminal 
wind velocity, and r' the position in the wind in units of stellar 
radius R*. From equation]^ we deduce a gas number density 
of ~ 1.7 X 10 ^ cm^'* at r = 8 R*, a value that is compatible with 
the number density of the pre-shock gas at 8 R* given in Table 
1^ Therefore, our modelled wind parameters are consistent with 
a wind acceleration fully complete at r ~ 8 R*. This radius cor¬ 
responds to the position at which the synthesis of silicate grains 
has ended. 

An increment of 10 in the gas number density pushes the ra¬ 
dius at which the wind is fully accelerated to r > 20 R*, accord¬ 
ing to equation and does not represent a model for the dust 
formation zone. However, such a number density enhancement 
may indicate trends produced by the presence of over-dense re¬ 
gions in the inner wind. 

In Figure 10 the grain size distributions for alumina at 
IR* and silicates at r = 6 R* are shown. We see that a denser gas 
produces much larger grains. For alumina, we consider the for¬ 
mation over one oscillation and no grain growth over several pul¬ 
sations. We see that most of the dust mass resides in large grains 
of size a ~ 0.2 /rm, while the size distribution peaks at a ~ 200 
A when no density enhancement is assumed (see Figure a. If 
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Fig. 10. Grain size distributions for a photospheric number den¬ 
sity xlO. Top; alumina in a shocked gas layer that has moved 
from 1 R* to 2 R*; Bottom: silicate in a shocked gas layer that 
has moved from 3.5 R* to 10 R*. 

deposition on grain surface were considered, grains could still 
grow between 0 = 0.5 - 1, to possibly start triggering a wind 
through Mie scattering very close to the star (IHofner 20081) . 

As for silicates, the grains reach an average size of ~ 0.1 
at r > 6 R*, with a dust-to-gas mass ratio of 5 x 10“^ at 6 R*. 
Because silicates form at larger radii than alumina, the impact 
of increasing the gas number density on the hnal grain size is 
weaker than for alumina, which forms in very dense layers at 
r < 2 R*. However, any inhomogeneity in the inner wind will 
result in the efficient synthesis of a population of large grains for 
both materials, and contribute to gas drag and wind acceleration 
close to the star. 

4. Summary and discussion 

Expanding the previous study by Duari et al. (1999), we present 
a comprehensive model of the inner wind of the M-type AGB 
star IK Tau, where the chemistry of the gas phase and dust clus¬ 
ters is coupled to the dust condensation phase, thereby providing 
a self-consistent treatment of the final dust masses and grain size 
distributions. The main results are as follows: 

- The modelled abundances for the molecules so far detected 
in the inner wind of IK Tau agree well with values derived 


from observations, except for PO and SO 2 (for the latter and 
NH 3 , the formation possibly occurs in a region located out¬ 
side the inner wind at r > 10 R*). While some molecules, al¬ 
ready present in the photosphere, are destroyed by the shocks 
and reform in the inner wind with abundances close to TE 
values, other species are only produced in the inner wind by 
shock-induced chemistry. The agreement for most molecu¬ 
lar species indicates that the periodic passage of shocks in 
the inner regions close to the star has a strong bearing on 
shaping the upper atmosphere and triggering the synthesis 
of molecules by non-equilibrium chemistry in AGB stars, as 
already suggested by Willacy & Cherchneff (1998), Duari et 
al. (1999) and Cherchneff (2016, 2011, 2012). 

- The prevalent gas-phase dust clusters in the nucleation phase 
are alumina and magnesium-rich silicates of forsterite sto¬ 
ichiometry (olivine). Enstatite (pyroxene) clusters are not 
abundant, with abundance values at least three orders of 
magnitude lower compared to forsterite. Small clusters of 
metal oxides (MgO and EeO) do form with low abundances. 
As for alumina, the clusters form very close to the star 
at r < 2 R*, while silicate clusters form at larger radii 
(r>4R*). 

- Eor our standard model of the IK Tau inner wind, one popu¬ 
lation of small alumina grains forms close to the star, peak¬ 
ing at a grain size of ~ 200 A. This population might be de¬ 
stroyed at each shock passage but will reform at a later phase 
in the pulsation. Therefore, our results indicate that a steady 
population of small alumina grains forms close to the star. 
The silicate population includes a component of large grains 
with a distribution that peaks at ~ 180 A after 6 R*. The 
mass of silicates is distributed in grains of smaller sizes be¬ 
cause of the lower gas number density at which the conden¬ 
sation takes place compared to alumina. When we consider 
an enhanced gas number density at the stellar photosphere, 
the size distributions of alumina and silicates are skewed to¬ 
wards larger grains, where alumina grains reach ~ 0.2 pm 
in size and silicate grains peak around 0.08 pm. Since we 
do not model the growth of grains through surface addition, 
these results can be taken as lower limits for the grain size 
and mass of alumina and silicates. 

- The total dust-to-gas mass ratios derived at 6 R* range from 
~ 1 X 10“^ and 6 x 10“^ for the standard and the enhanced 
density cases, respectively. This is in good agreement with 
dust-to-gas mass ratios derived from observations for IK Tau. 


Several recent interferometric studies of oxygen-rich AGB 
stars point to the presence of an alumina shell close to the star. 
While Norris et al. (2012) found a halo of large (~ 0.3 pm), 
transparent grains close to the photosphere (at r < 2 R*) in a few 
O-rich AGBs, Zhao-Geisler et al. (2012) reported evidence for 
the coexistence of a warm layer of H 2 O and a AI 2 O 3 dust shell 
close to the stellar surface in several O-rich evolved stars. The 
presence of two dust shells of alumina and silicate composition, 
where the alumina shell is located below the silicate shell, have 
been detected in a few O-rich stars (IZhao-Geisler et al. 20121 
IKarovicova et al. 2013b . This latter interferometric study points 
to the formation of alumina dust at ~ 2 R* and of silicates 
at ~ 5 R* from the stellar surface for the few sources show¬ 
ing the typical alumina and silicate features in the mid-infrared. 
Therefore, these various studies point to a dust condensation sce¬ 
nario similar to that derived from the present study, with possibly 
two separate dust populations, one close to the star and consist¬ 
ing of alumina grains, the other further out and made of silicates. 
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Owing to the small opacities of iron-free silicates in the 
wavelength range where most of the stellar emission flux is pro¬ 
duced (around 1 fim), acceleration of the gas by photon scat¬ 
tering on large forsterite grains was proposed as a substitution 
to photon absorption (IHofner 20081 IBladh et al. 20121 2013). 
According to these studies, Mie scattering on grains in the size 
range 0 . 1-1 is able to generate a wind acceleration and 
terminal velocity, and mass loss rates that agree with values de¬ 
rived from observations. Whether such large forsterite grains re¬ 
ally form in the wind requires a careful chemical kinetic-based 
investigation, because models by Hofner (2008) and Bladh et 
al. (2012, 2013) do not treat the nucleation of forsterite grains 
but describe growth by assuming a pre-existing populations of 
small seeds in the inner wind. The interferometric study of a 
few oxygen-rich AGBs by Norris et al. (2012) indicates scat¬ 
tered light dust shells very close to the star, at radii < 2 R*. An 
iron-free silicate grain composition is favoured because iron-rich 
silicates would be too hot to survive so close to the star owing to 
their high near-IR opacity, but alumina grains are not excluded 
and can equally well reproduce the observations. The derived 
average grain sizes are large for all objects, and around 0.3 ^m. 
Our chemistry-based model of IK Tau shows that silicates of 
any type should not form at radius < 3.5 R*, because the for¬ 
mation of small molecular clusters involved in the nucleation of 
silicates is hampered. On the other hand, alumina dust clusters 
can easily form close to the star from simple three-body associ¬ 
ation of AlO molecules, and further oxidation to form the small 
ring AI 2 O 3 , and the subsequent alumina dimers at gas tempera¬ 
tures less than 2500 K. The high gas density warrants a rapid and 
efficient condensation of these alumina clusters to form grains, 
whose total mass is distributed in a population of medium-size 
grains with radius ~ 0.02 jjm. Possible surface growth will take 
place when the post-shock temperature is low enough to allow 
the grain survival and growth, which will result in larger grains. 
These results and trends then support the presence of an alumina 
shell very close to the stellar surface, as observed by Norris et 
al. ( 2012 ). 

The formation of silicates takes place at larger radii (r > 
3.5 R*) because the nucleation scheme involves the formation of 
several molecular intermediates, specihcally the silaformyl radi¬ 
cal, HSiO, and silanone, H 2 SiO, whose dimers cannot form close 
to the star. Our model for IK tau gives silicate grains of medium 
size in the range 200 A - 0.1 jum, and most of our forsterite 
grains are formed between 6 R* and 8 R*, at which point con¬ 
densation stops because the local gas density is too low to facil¬ 
itate efficient coalescence and coagulation. In this radius range, 
the wind will be accelerated and this process will gradually 
freeze the gas-phase chemistry and the formation of molecules 
and dust clusters. It is thus unlikely that the dust shells de¬ 
rived from interferometric studies, and found close to the stellar 
surface, consist of iron-free silicates (e.g., ISacuto et al. 2013t . 
In our model, grains of alumina are formed close to the star. 
Because they are likely destroyed at each shock passage and re¬ 
form in the post-shock gas, they will reside close to their forma¬ 
tion location without being expelled at larger radius. If so, and as 
mentioned before, our results indicate the presence of two dis¬ 
tinct grain populations in the inner wind, as recently detected 
by Karovicova et al. (2013). We note that the various interfero¬ 
metric observations often relate to semi-regular objects and not 
variable Mira-type stars like IK Tau. Because semi-regular AGB 
stars usually have large photospheric temperatures, the forma¬ 
tion of silicates may be delayed to even larger radii, and thus 
much smaller amounts of small silicate grains will form under 
more diffuse gas conditions. Therefore, alumina might be the 


major dust component in semi-regular variables whereas silicate 
dust might become the prevalent dust type in the cooler and more 
evolved oxygen-rich AGB stars. 


As already mentioned, we do not consider the growth of 
dust grains through surface addition of gas-phase species. It is 
often assumed that some small seeds formed close to the stel¬ 
lar surface could act as condensation nuclei and grow from 
heterogeneous addition of material on the surface. Specihcally, 
Ti 02 has often been proposed as a condensation seed in AGB 
winds and brown dwarf clouds because of its refractory na¬ 
ture (Jeong et al. 2003 ILee et al. 2015t . Furthermore, Goumans 
& Bromley (2013) showed that one Ti 02 molecule was enough 
to kick-start silicate dust formation via SiTiOa seeds, thus over¬ 
coming the SiO dimerisation bottleneck. Our results on gas- 
phase TiO/Ti 02 indicate that titanium remains in atomic form 
and in TiO in the inner wind, with some TiO converted into 
Ti 02 by shock chemistry. Several transitions of molecular TiO 
and Ti 02 have been observed in the wind of the supergiant 
VY CMa (IKaminski et al. 2013al IDe Beck et al. 2015]) . where 
the Ti 02 emission is extended and not always correlated with 
the dust continuum emission. These results point out that Ti 02 
clusters and grains may not play an important role as condensa¬ 
tion seeds. It seems that the modest amount of Ti 02 we derive in 
the gas phase precludes solid Ti 02 to be a key agent in the dust 
formation process in IK Tau and other O-rich AGBs. This con¬ 
clusion was also drawn by Plane (2013) in his study of calcium 
titanate formation in AGB stars. 


The formation of iron-rich silicates is not considered in the 
present study because there are currently no established chem¬ 
ical nucleation routes that lead to the formation of such a ma¬ 
terial at high temperatures. The iron and magnesium abun¬ 
dances are very similar in the stellar photosphere, and iron 
seems not to be depleted in other dust types in the inner wind. 
Indeed, we do not form the molecules FeO and Fe 2 in signif¬ 
icant amounts, two molecules that are essential to trigger the 
formation of iron oxides and pure iron grains. Most of the iron 
is in atomic form for r < 10 R* and might then easily be in¬ 
cluded in small silica-based clusters, as is the case for magne¬ 
sium, to form Fe-rich silicates. If so, two silicate-based compo¬ 
nents are expected, one Mg-rich and the other Fe-rich, as was ob¬ 
served in laboratory study of a mixed Mg-Fe-Si 0 -H 2-02 vapour 
( |Rietmeijer et al. 1999) . Because iron-rich silicate grains have 
high opacities around 1 //m, they will be more efficient at driving 
a wind through radiation pressure and hotter than their iron-free 
counterparts, and will thus favour the emission of typical silicate 
features at 9.8 i^m and 18 fim usually observed in the SEDs of 
Mira-type stars. 

The inner wind gas conditions may be conducive to the syn¬ 
thesis of other solid compounds like spinel, MgAl 204 , and sil¬ 
ica, Si 02 . Because of the efficient condensation of alumina very 
close to the star, and the availability of atomic magnesium (MgO 
does not form in large amounts at r < 2 R*), the formation of 
spinel may occur along with that of alumina, or possibly com¬ 
pete with the latter. Furthermore, we have already mentioned 
that the production of alumina is triggered by the dimerisation of 
AlO. However, the large abundance of AlOH in the inner wind 
may induce the formation of molecular species not considered in 
the present study, such as HAIOH. Enhanced AI 2 O 3 production 
may occur through HAIOH dimerisation and processes mimick¬ 
ing the proposed routes through silanone considered for silicate 
formation. Such chemical pathways have not yet been studied 
theoretically or experimentally, and we plan to investigate these 
new chemical routes and their impact on alumina production. 
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combined with the formation of spinel and iron-rich silicates, in 
a future study. 
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